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, Abstract 

The goal of this paper is to exhibit a critical mass phenomenon occuring in a model for cell self-organization via chemotaxis. 
The very well known dichotomy arising in the behavior of the macroscopic Keller-Segel system is derived at the kinetic level, 
being closer to microscopic features. Indeed, under the assumption of spherical symmetry, we prove that solutions with initial 
, data of large mass blow-up in finite time, whereas solutions with initial data of small mass do not. Blow-up is the consequence 
of a virial identity and the existence part is derived from a comparison argument. Spherical symmetry is crucial within the 
two approaches. We also briefly investigate the drift-diffusion limit of such a kinetic model. We recover partially at the limit 
the Keller-Segel criterion for blow-up, thus arguing in favour of a global link between the two models. 
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+3.! Introduction 

a . 

Qn| this paper we aim to exhibit a blow-up versus global existence phenomenon for a kinetic model describing collective motion 
1 oi cells in two dimensions of space. The so-called Othmer-Dunbar-Alt system [2"5l |3"T] reads as follows, 

dtf + v-V x f=[ T[S](t,x,v,v')f(t,x,v')dv' -X[S](t,x,v)f(t,x,v), x&R 2 , t>0, 

Jv ' eV (1.1) 



— AS + aS — p(t, x) — j f(t,x,v)dv, 



L_w|iere f(t,x,v) denotes the cellular density in position x velocity space, and S(t,x) is the concentration of the chemoattractant. 
^j3Tie velocity set is assumed to be bounded, and for simplicity we take V — B(0, R) throughout this paper. The initial data fo 
(^j5>eiongs to i 1 (R 2 x V) (more appropriate assumptions on / will be stated within Assumption Theorems [2] and [4]) . Observe 
^that the mass of cells is formally conserved in time: 

•l-H 

X 



f(t,x,v)dvdx= // fo(x,v)dvdx. 

xV J JR 2 x V 



The turning kernel T[S](t, x, v, v') > denotes the probability of transition between velocities v' — > v at position x and time 
t, and X[S] — J v , eV T[S](t, x, v', v) dv' is the intensity of this Poisson process. The influence of the chemical field S is highlighted 
in the notation T[S}. We make the particular choice 

T[S](t,x,v,v')= X o{v-VS(t,x)) + , (1.2) 

with constant xo > 0, that is to say cells choose only favorable directions when they reorientate, and they align more likely 
with the gradient of the chemical. This mechanism is actually not well-suited for describing accurately bacterial motion like 
the 'run and tumble' process performed by E. coli [13) . However, this fits well with motion of bigger and more complex cells 
capable of sensing a space gradient of the chemical and to orientate accordingly (amoeboid or mesenchymal motion |14j). More 
complex kinetic models involving saturation effects or interactions between cells and the surrounding tissue have been investigated 
respectively in [TU1 [H] . 
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1.1 Statement of the main results (blow-up vs. global existence) 

The kinetic model under consideration in this work can be written more precisely as follows: 

| dtf + v ■ V x f = Xo(v ■ VS)+p- X oLo\VS\f, v€V = B(0,R), i£R 2 , t>0, 
\ -AS + aS = p[t,x). 

where uj = J v , eV (v' ■ V5/| V5|) + dv' = 2i? 3 /3 (see Lemma [HI below, the precise value of u> does not play an important role in 
the sequel). 

Assumption Al (Initial datum). Assume that the initial density fo{x, v) > belongs to L 1 f] L°°(R 2 x V). Assume in addition 
that fo has spherical symmetry: for any rotation ®, fo(@x,@v) = fo(x,v). We denote by M the total mass: 

M = // fo(x,v)dvdx. 

J JR 2 x V 

This assumption ensures that the system (|1.3[) has a unique solution, which remains spherically symmetric during its life 
span (see Section [2] and Appendix). 

Definition 1 (Blow-up in the context of kinetic chemotaxis). A solution of (|1.3|) is said to blow-up if after some time (possibly 
infinite), f(t, x, v) exits LV X L1 for all exponents p and q such that 

2<p, Kq, o<--i<i. 

q p 2 

This particular choice of exponents is clearly related to the dispersion Lemma\^ In fact L V X L\ turns out to be a natural space for 
existence theory associated to (|1.3[) (see Appendix and fljjj). 

Theorem 2 (Blow-up for large mass under spherical symmetry, a — 0). Assume that fo satisfies Assumption and has a 
finite second moment w.r.t. the space variable x. Assume that the mass is large enough: 

Then the solution of (|1.3[) with a = blows-up in finite time. 

Corollary 3 (Blow-up for large mass under spherical symmetry, a > 0). Assume that fo satisfies Assumption J$l\ and has a 
finite second moment w.r.t. the space variable x. Assume that the mass is large enough (jl.4|) . Assume in addition that the 
following condition is fulfilled initially, 

\x\ 2 p (x)dx<C(xo,M 1 \V\), (1.5) 

R 2 xV 

where C{xotM, \V\) is an explicit (but heavy) constant vanishing when the case of equality is reached in (jl.4p . Then the solution 
of (|1.3|) with a > blows-up in finite time. 

Before we state our next result we introduce some notation. For any exponent < 7 < 1 define: 

1 r /2 

0(7) = l + ~ / |sin6>n 7 d6. (1.6) 

7T Je = -7T/2 

Theorem 4 (Global existence for small mass under spherical symmetry). Assume that fo satisfies Assumption J$Ji and lies below 
fco|a;| _7 for some < 7 < 1 and ko > 0. Assume in addition that the mass is small enough: 

M < ^1 . (1.7) 

Then the solution of (|1.3[) exists globally in time. 

In fact we can derive some pointwise estimate, and get that the solution does not blow-up in infinite time either. 

Corollary 5 (Simplified criterion for global existence). Assume that fo satisfies Assumption ^[7] and additionally that fo € 
C"(K 2 ). There exists a threshold M*(xo, \V\) such that if M < M* then the solution of (|1.3|) exists globally in time. 
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The two theorems above and their corollaries require a few comments. 

(i) The meaning of spherically symmetric solutions in the context of kinetic models is given in the Appendix. 

(ii) The case of the velocity space being V = §(0, R) is also investigated as a direct adaptation of Theorem [2] (see page[TT|). We 
obtain that blow-up occurs if the criterion (|1.4[) is replaced with: M > (16n/xo\V\). In particular, cells' velocities are not 
likely to converge to zero in the general case when blow-up occurs (both density and velocity collapse) and blow-up can 
happen even if the set of admissible velocities is bounded from below. In this situation, the tumbling frequency dramatically 
increases in the neighbourhood of the blow-up location. 

(hi) Note that Theorem^s included in its Corollary [3] when a — > 0. However, we state first Theorem[3]for the sake of coherence 
(the Corollary can be thought of being a perturbation result for a > 0). 

(iv) We shall derive in Section [3] a weaker criterion for blow-up when a > 0, involving the initial first and second moments, and 
the initial current (|3.13[) . However we opted for a simplified presentation in Corollary [3] involving the second moment only. 

(v) Theorem 2] and its corollary [5] are concerned with a = 0, but the case a > can be proven similarly. 

(vi) In Theorem [3] the assumption fo(x,v) < fcoM -7 is not very restrictive, because the tail \x\~~ 1 is not integrable at infinity. 

(vii) We might have expected < 7 < 2 instead of < 7 < 1 because we are in dimension 2, and the natural space for solutions 
here is L 1 (R 2 x V"). However it appears in (|1.6[) that 7 < 1 is a crucial condition that we are not able to overcome (it 
underlies the fact that our reference function k(x,v) we are comparing with in Section [4] is not integrable with respect to 
velocity if 7 > 1). 

(viii) The continuous function 7/^(7) goes to zero both for 7 — > and 7 — ► 1. Therefore there is a best compromise 7* which 
maximises the condition (jTTTJ) : M < M* = 27*/(xo|P r |fi(7*)). However in this paper we keep general 7 for the sake of 
clarity. Corollary O is nothing but applying Theorem U to 7 = 7* defined as above. 

(ix) Notice that the two critical mass thresholds (resp. (|1.4[) and M*) do not match, because fl > 2 ensures M* < / {xo\V\) < 
27r/(Xo|V|). Numerically we obtain 47*/^( 7 *) « 0.806 < 32! 

To the best of our knowledge, few blow-up results have been exhibited for kinetic models. Let us mention primarily the virial 
identities derived by Horst [55] , Glassey and Schaeffer [TH] respectively for the Vlasov-Poisson and the relativistic Vlasov-Poisson 
models in the gravitational (i.e self-attracting) case (see also [18] for a presentation of these results). These remarkable identities 
allow the authors to show blow-up of the solutions having negative energy (resp. in dimension d > 4, and in dimension d = 3 
under spherical symmetry) . Recent progress aims to describe precisely the blow-up dynamics using the Hamiltonian structure and 
concentration compactness techniques for the Vlasov-Poisson system [33]. Within the context of chemotaxis models, Chavanis 
and Sire have derived various virial theorems which share several common features with the identities derived in this paper [TT] . 

Singularity formation plays a very important role in the kinetic theory of Bose-Einstein condensates, which arise when part 
of the particle density concentrates in the same quantum state jT51 [3H [25] • No virial identity has been developed in this theory 
however (in fact it would be irrelevant to argue through a vanishing second moment for proving concentration in this context): 
the convergence of the solution towards a singular limit (a Dirac mass together with a regular part) at low temperature for 
the Boltzmann-Nordheim (resp. Boltzmann-Compton) equation is performed via entropy/entropy dissipation techniques (see 
p~5l [26] and [32] for an overview of the kinetic theory of Bose-Einstein condensates) . Virial identities in the context of kinetic 
theory are also developed in cooling processes within the Boltzmann equation where particles are subject to inelastic collisions 
[3J[37l[3]. However, the context of the last two examples differs from the situation we are interested in because concentration 
occurs in the 'velocity' variable whereas for cell chemotaxis and self-attracting Vlasov-Poisson systems it occurs in the space 
variable. Consistently enough, the two aforementioned examples deal with homogeneous in space kinetic equations. 

The paper is organized as follows: the next introductory subsections enlarge the picture, and replace the above results in a 
more general framework. We show the strong continuity between this critical mass phenomenon and previous existence theorems 
in kinetic theory for chemotaxis, and we explain the strong links between the current kinetic model, and the so-called parabolic 
Keller-Segel model. In Section [2] we briefly state a local existence theorem, which is later given with full details in Appendix. 
Section [3] is devoted to the proof of the blow-up result having different variants, using an efficient virial identity. In Section [4] we 
prove global existence by a comparison argument with a very specific (and singular) reference function. The Appendix contains 
a short description of the meaning of spherically symmetric solutions in the context of kinetic equations, and also a general local 
existence Proposition (based on dispersive estimates) to be crucially used in the comparison argument of Section [U 
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1.2 Brief review of context in the light of existence theory 

As a by-product of this critical mass phenomenon we can argue that previous existence results in the kinetic theory of chemotaxis- 
biased cell motion were not far from being critical. The state of the art of existence results is discussed in the introduction of 
5 . We recall below some of the issues for the sake of puting the results of our paper in context. From (jl.lj) two classes of 
problems emerge depending on the assumption conditioning the turning kernel. The same tool comes out to be powerful in both 
situations. As it will be used at few places in the present paper, it is worth recalling the dispersion lemma [7] which measures 
the action of the free transport operator on mild L^L^-norms: 

Lemma 6. (Dispersion estimate) Let go(x, v) 6 L q (M. 2 ; L P (M. 2 )) where 1 < q < p < oo, and let g solve the free transport equation 

d t g + v-V x g = 0, (1.8) 

with initial data g(0,x,v) ~ go(x,v). Then 

y^hgLt < f2 (i/g-i/p) iisoHlsj* ■ ( L9 ) 

Strichartz estimates [7] have also been shown to apply succesfully to those run-and-tumble problems (see [4] and Remark [8] 
below) . 

Transport-dominating regime. It deals with the case where the turning kernel can be estimated pointwise in terms of 
Sobolev norms of the chemical signal. For instance assumptions 

T[S\(t,x,v,v') < C\\VS(t)\\ 1 - u , < v < 1, or T[S](t, x, v, v') < C\\VS(t)\\ r , r < oo , (1.10) 

both lead to global existence of solutions (some superlinear power can in fact be added in the second case, see [5] for details). 
When estimating the evolution of LgL^-norms, the dispersion due to the free transport operator turns out to have a strong 
enough effect to counterbalance the aggregation due to the tumble kernel. 

Derealization effects. It deals with the case where the turning kernel is pointwise estimated through some space dereal- 
ization (volume effects, protrusion sending), for example, 

T[S](t,x,v,v') <C\VS(t,x + ev)\, e > 0. (1.11) 

In this case the solution is again proven to be global in time [23]. In [4] a second derivative subject to the same sort of 
derealization can even be added to (|1.11|) . 

It is worth noticing that our special choice of turning kernel (|1.2p is critical for both assumptions p.!0[) and p. lip . 

Remark 7 (Dispersion method is borderline.). The turning kernel we are studying in this paper satisfies 

T[S](t, x, v, v') < C ||VS(t)|| £00 • (1.12) 
It is natural to ask whether this property alone implies global existence. Indeed, if we work as in '4, 8j we arrive at 

\\pit)\\LlL% < J \\VS{s)\\ LX \\p{s,x-{t-s)v)\\ LlLi da. 
where 1 < q < p < oo. By Lemma® the right-hand side can be controlled by 

C(y) I s 2(i/1-i/p) I' VS ^h- \\P(S)L« ds ■ 
For p we can use interpolation: \\p\\ L g < ||p||^i P ^ 9 IIpIIlp 7 ' 9 • For V S we can use the elliptic estimate (see Lemma \7l\ ): 

Wvsw^kcWpWI/^WpwH 2 ,2< P <oo. 

We obtain eventually 

\\p(t)\\iZLZ<cf* g2(1/ ^ 1/p) \\p(t-s)\\il 2+p ' /9 ' ds, 

and it is impossible to achieve both 2 (1/q — 1/p) < 1 for integrability near s = and p' /2 + p' / q' < 1 for applying a global 
Gronwall's lemma. 



4 



Remark 8 (Strichartz method is borderline.). If we are willing to impose a smallness condition we can try to use Strichartz 
estimates in the same spirit as [J^l. Recall from J^j that if f solves the free transport equation 

dtf + V- X7 x f = g, 

then 

Wf\\LiL% Ll < Co + CiNI Lf ' isLS • C 1 - 13 ) 

where Cn depends only on the initial data and the parameters q, p, r satisfy 

l<r <P<oo, - = 2(---)<l, - + ->1. (1.14) 
q \r p J p r 

The borderline case (when the middle condition in (|1.14| fails) would correspond to q = 2, p = 2, r = 1, 

WfW^LlLl < CO + Cl hh^Ll ■ (1-15) 

We are interested in g = C \\VS\\ Lao p under the assumption of spherical symmetry. In this special case we have 

1/2 / r r \ 1/2 



Vx \VS(x)\ = -_l p(X)XdX<^U XdXj p{X) 2 XdX) <C||p|| z , 



i r i 

- / p(X)X dX<- 
r .In r 

therefore 

Wf\y t LiLi< Co + Ci\\\\p(t)\\L* pitrfW^^Co + CtM \\f\\ LlLlLl . (1.16) 

// the mass M was small enough we would be able to bootstrap. 

Thus an alternative proof of global existence under small mass and spherical symmetry (much simpler than the one we develop 
in Section^ would rely on a critical Strichartz estimate that we are not currently able to handle. 



1.3 A reminder of the classical Keller-Segel in 2D of space 

The critical mass phenomenon studied in this paper shares several similarities with the qualitative behaviour of the parabolic 
Keller-Segel system in two dimensions of space: 

f d t p = Ap- xoV • (pVS) , t > , x e R 2 , 

{ (1-17) 
[ -AS + aS = p. 

In fact, there is a simple dichotomy: if the mass is below the threshold M < 8n/x then the solution is global in time and disperses 
with the space/time scaling of the linear heat equation; on the other hand, if it is above the same threshold M > 8ir/x, then 
the solution blows-up in finite time (in the case a = 0). For blow-up in the case a > one usually adds an hypothesis close to 
(|1.5[) [6] . This critical mass phenomenon was first derived in a bounded domain with radial symmetry [24l [28] . Energy methods 
based on ad-hoc functional inequalities (either Trudinger-Moser or Hardy-Littlewood-Sobolev with a logarithmic kernel) were 
developed later on [T71 Q] . 

The analogy is not complete however, as can be seen in the details of our proofs. Concerning the blow-up, we have to 
differentiate twice in time the virial identity as opposed to Keller-Segel for which it holds true (when a = 0): 



d 1 
dt 2 



R 2 



\x\ 2 p(t,x)dx = 2M (l-^J 



Concerning global existence, the Keller-Segel system is equipped with a free energy (entropy minus chemical potential energy) 
which is dissipated along the trajectories and this provides useful a priori estimates ensuring global existence for small mass. 
No such energy is known at the kinetic level. In the present work we use a comparison principle with a singular but integrable 
reference function. 



1.4 Drift-diffusion limit (formal) 

The parabolic Keller-Segel system can be obtained as a drift-diffusion limit of the kinetic Othmer-Dunbar-Alt model [30l [9l [8]. 
when the chemotaxis bias is a small perturbation of an unbiased process. We may express this fact by modifying the turning 
kernel under consideration: 

T e [S](t,x, v, v') = F(v) + e X o(v -VS(t,x)) + , (1.18) 

instead of (|1.2[) . Thus, the jump process consists in the superposition of a relaxation process (towards a velocity distribution 
F(v) > such that J vF(v)dv = and J F(v)dv — 1) and a small bias due to chemotaxis. We assume F(v) to be rotationally 
symmetric (in order to match with the context of this paper). 
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Remark 9. Previous works (see e.g. 176) /] state in general that the scattering operator T e [f,S](x,v) — J v T £ [S]f(v')dv' — 
X e [S]f(v) can be decomposed as To[f] +eTi[f, S], where the linear unbiased operator T~o[f] possesses an equilibrium configuration 
with respect to velocity, namely there exists a probability distribution F(v) such that J v vF(v)dv = and Tq[F] = 0. In this 
paper however we restrict the presentation to the special case of relaxation towards F for the sake of clarity, but the standard 
procedure can be performed in the same way. 

The kinetic model with the parabolic scaling writes 

ed t f e +v-V x f e = - ( Pe F(v) -f e + e Xo (v ■ VS e ) + p £ - e X0 Lo\VS £ \f e ) , 

^ (1.19) 

— AS e + aS e = p e (t, x) . 

Formally, as e — ► 0, the cell density f e decouples into a product p(t,x)F(v) (so that the leading order term cancels), and p is to 
be determined. To do so, integrate against 1 and v the first line of (|1.19[) and get respectively 

edtPe + V • j e = , with j e = / vf e dv , 

Jv 

edtj e + V • (J v <g> vU dv^J = ^ (-j e + £Xo^ (VS e ) Pe ~ e X ouj\VS t \j^ . 
Still formally, we obtain the renormalized flux for small e 

j = -V • QT v ® vf e dv^j + Xo ^ (V5 e ) p t . 
Therefore we obtain as e goes to zero, the limiting parabolic equation for the cell density in space p(t, x), 



d tP = V 



v (g> vF(v) dv 



v 



V " ) - ' (^) , (1-20) 



coupled with the chemical potential equation —AS* + aS = p. 

Notice that the current assumptions fit with the framework of [S], which makes this analysis rigorous for short time t < t* 
(independent of e) . We refer to the end of Section 13.31 for a discussion about this formal limit from the viewpoint of blow-up 
results. In a short, we show that blow-up criterions are indeed the same (asymptotically) for the kinetic model (| 1 . 3[) and its 
parabolic limit. This raises the question whether this convergence is still valid for larger times. In other words: do the kinetic 
and the parabolic Keller-Segel systems remain close to each other throughout their respective periods of existence? Our result 
strongly supports the fact that the solutions are indeed close for all time (before the blow-up time). However, a rigorous statement 
together with a full proof of convergence for all time has yet to be performed. 



2 Preliminaries: local in time existence and uniqueness 

In this Section we prove local existence and uniqueness (without the assumption of spherical symmetry) for the system (|1.3p . 

Proposition 10. Consider the model (| 1 - 3[) with the turning kernel given by (II. 2\ . Fix p S (2,oo) and suppose that /o € 
h\ v n LP v . Then there exists a positive number T depending only on Jo an d a unique solution f with 

f el°° ([0, T]-L l (R 2 x V)) n L°° ([0, T];L P (R 2 x V)) 

Before going into the proof of Proposition [Hjl let us state a useful elliptic estimate. We omit the easy proof which is a direct 
consequence of the Hardy-Littlewood-Sobolev inequality. 

Lemma 11 (Elliptic estimate). Let p > 2 and S be the solution of —AS" = p in the sense 

2tt J \x-y\ z 

Then 

HVSHoo < C( P ) Wpf-Z \\p\\J P , lim C(p) = +co . 



G 



Note that this elliptic estimate holds true for p = 2 in the spherically symmetric case. However, we shall not use that variant 
here. 

Proof of Proposition \1(A Let us write the nonlinear kinetic equation of interest as 

dtf + v-V x f = N{f), 

where the nonlinear scattering operator is given by 

N(f) = X o («-V5) + p-xow|VS|/. 
We shall prove that the nonlinear operator satisfies a Lipschitz estimate, 

\\N(h) - N(f 2 )\\ x < L , ||/ 2 || x ) ||/! - / 2 || x , (2.1) 

where the norm is defined by 

h\\x = sup (\\g(t,x,v)\\ L1 + ||ff(t,ar,«)|| iS , I ) . (2.2) 

0<t<T v x '" ' ' 

We split the difference of the nonlinear contributions into four different parts, namely, 

N(fi)-N(h)= Xo((vVS 1 ) + -(vVS 2 ) + )pi (I) 
+Xo (v ■ V5 2 ) + ( Pl - p 2 ) (II) 
-Xow(|V5i|-|V5 a |)/i (HI) 
-X0WIVS2KA-/2) (IV). 

For the first contribution I we have 

\I(t,x,v)\ <C( X o,V)\\VS 1 (t)-VS 2 (t)\\ LOO \ Pl (t,x)\. 
The difference of the two gradients in L°° can be estimated via the elliptic estimate of Lemma [TTJ 

||VS!(t) - VS 2 (t)\\ LX < C( P ) \\pi(t) - P2 (t)\\ 1 L^ \\pi{t) - P 2{t)\\i , 



<c(p,v)( sup ||/i(i')-/ 2 (Oil li 2 su p Wh(t')-f2(t')\\i P 

\0<t'<T J \0<t'<T 

<C(p,V)\\h-h\\ x . 

As a consequence we get a Lipschitz condition for the first part /: 

|I(i,a;,tO| < C(p,xo,V)||/i-/ 2 || Jf |pi(t,o;)|, 
\\I\\ Llv < C(p,x a ,V)\\h~h\\ x \\ Pl \\ LP 

< C(p, X o,V)\\fi-f2\\x\\fi\\x- 

Similarly 

||/|| ii>B <C(p,Xo j VO||/i-/ 2 || a -||/iL . 

therefore 

imix<C(p,Xo,y)||/i|| x ||/i-/ 2 || x . (2.3) 

The estimates for 77, /// and IV are obtained analogously aso that we end-up with the desired estimate (|2.f p with L (||/i || x , ||/ 2 || x ) = 
C(j>,xo,V)Q\fi\\ x + \\M\ x ). 

To conclude, let us mention that the norm || • \\x defined by (|2.2p is preserved through the action of the free transport 
operator, thus a fixed-point argument can be developped and leads to the conclusion. □ 

Remark 12. In the appendix we state a more complex existence/uniqueness result in suitable spaces L V X L\. This is to fit with 
the comparison method of Section^ 
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3 Formation of a singularity for large mass 

3.1 A blow-up criterion in the case a = 

We first need to state a technical Lemma for explicit computations. 
Lemma 13 (Averaged quantities). Recall that V = B(0,R). 
(i) For any g£K 2 we have 

J (v ■ q)+dv = ^j-M ■ 

(ii) For any (p,q) £ K 2 x I 2 , we have 
Proof. Item (P is immediate. 

Concerning item (JTTJ) , denote J(p, q) — J v€V {p ■ v)(v ■ q) + dv, then J is symmetric: 



J(p,q) = I (p'v)+(vq) + dv- (p'v)-(v-q)+dv 

(p ■ v) + (v ■ q) + dv — / (—p-w)-(—w-q)+dw 

J J w 

(p ■ v) + (v ■ q) dv . 

Moreover, J(j>,q) is linear w.r.t. p, so it is bilinear w.r.t. (p, <?). It remains to compute the associated quadratic form: 

J{P,P) = / (p-v)(v -p)+dv 
= 2 ( {P ■ v){v ■ p) dv 



Thanks to isotropy we obtain, 



R /-27T 

2 „„„2 , 



v®vdv = II / r z cos? 6 rdrdO \ Id (3.1) 

= 7T— — Id . 

4 



Consequently we deduce 



J(p,q) = —^-(p ■ <?) • 



□ 



Proof of Theorem^ We plan to evaluate explicitly the time evolution of the second moment w.r.t. to space variable x. We 
introduce the notation, 

I ( t ) = o / / M 2 /(*> x i v ) dvdx > 

l 2 xV 



2 
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We differentiate twice in time 
d 



jJtt) = / (x-v)f(t,x,v)dvdx 
1 

+ 2 



x \ 2 \j J T[S](t,x,v,v')f(t,x,v')dv'dv- I \ T[S](t,x,v',v)f(t,x,v)dv'dv\dx 
(x ■ v)f(t, x, v) dvdx , 



' x J V 

|2 

\V 

t X J V 



jjal{t) = I I \v\ 2 f(t,x lV )dxdv 



(x-v)T[S](t,x,v,v')f(t,x,v')dv'dvdx — / / (x ■ v)X[S](t,x,v)f(t 7 x,v) dvdx 



x J v J V' 



With the particular choice for T[S] given by (|1.2p (it does not depend on the anterior velocity v'), we obtain: 
^jJ(t) = JJjv\ 2 f(t,x,v)dvdx 

+Xo / / {x ■ v)(v ■ \/S)+p(t, x) dvdx — xo / I i x ■ v) < I {v' •'VS)+dv'> f(t,x,v)dvdx 



Therefore, applying Lemma [T3l we get 

^I(t)= [ [ \v\ 2 f{t,x,v)dxdv + X o^- fx- VS{t,x)p(t,x)dx- X o^- f [ (x ■ v)\ VS|(t, x)f(t, x, v) dvdx . (3.2) 

The following computation is well-known within the theory of the Keller-Segel system pQ : 

1 f f x — y 
x ■ V S(t,x)p(t, x) dx = — — I I x ■ -. prp(t,y)p(t,x) dydx 



2tt J x J y \x~y\< 
1 f ft \ X ~V 



Therefore we obtain from (13.2 



i i i (x-y)- | ToP{t, y)p(t, x) dydx 

^JxJy \x-yr 

47T 



^I(t) < i?2 M -^M 2 - X0 ^ J J(x-v)\VS\(t,x)f(t,x,v)dvdx 

< ^ M (l-^^^-XD^JjB-j(t,x)\VS\{t,x)dx, (3.3) 
where the current j(t, x) = J vf(t, x, v) dv satisfies 

dtp + V • j = . 

Introduce the notation, 

j = fi > M (*2|^-l), (3.4) 

which is positive by assumption (| 1 .4[) of Theorem [2j 

In spherical coordinates, we can compute exactly the contribution of the remaining term in (13 . 3|) , using the following identities 
(see Appendix): 

\VS\(t,x) = \S'(t,r)\, r\S\t,r)\= f Xp(t,X)dX, 

J\=o 

i(t,x) = j%r)-+^(i,r) — , r j\\(t,r) = -— Xp(t,X)dX=— Xp(t, X) dX . 
r r dtj x=0 OtJ x=r 
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Therefore, under the hypothesis of spherical symmetry we get 

x ■ j(t,x)\WS\(t,x)dx = -2vr / rj 11 (t, r)\S'(t, r)\ rdr 

Jr=0 



X=r 



d r°° r 00 d r°° ( r°° \ 2 

-M— / Xp(t,X)dXdr + ir— / Xp(t,X)dX) dr 

dt J r= o J \ =r dt J r= o \Jx =r J 



Integrating once in time the inequality (|3.3p leads to 



where K(t) is defined by 



d d 27? 3 

- m <^J(t)\ -H + xo— (*(0) -*(*)) 



OO pOD pOO / pOO \ 2 



K{t) — Ml / Xp(t, A) dAdr - tt / / Ap(i, A) dX dr 

Jr=0 J\=r Jr=0 \J \=r J 

Xp(t, A) dXdr + tt f — - / Xp(t, X)dX) ( / Ap(i, A) dA ) dr . 



r=0 



Thus if(i) is obviously a nonnegative quantity. It is worth noticing that K (t) is finite provided that the density p has a finite 
second moment. Indeed we have 



p OO />oo 

K(t) < M / Xp{t,X) 

Jr=0 J \=r 



oo />oo 

dXdr 



OO 

2 , 



= M / X z p(t,X)dX 
J\=o 
M f 

< 7T / \x\p{t,x)dx. 



We deduce in particular using the Cauchy-Schwarz inequality, 

K(t) < ^-M 3 / 2 ^/2I(J). (3.5) 

Remark 14. It seems surprising that the contribution of the loss term to the virial identity can be written as the derivative 
of a nonpositive quantity. It is not surprising however if we notice that in spherical coordinates, this contribution is the scalar 
product between rj"(t,r) and rS' (t,r) satisfying respectively 

d 

On the other hand this property still holds true when the turning kernel contains a linear part (see (|3.14[) ). 
We end up eventually with 

d r f 2R 3 

T I{t)< (x-v)fo(x,v)dvdx-St + X o^ r K(0). (3.6) 



This proves that the second moment formally vanishes in finite time, therefore a singularity necessarily forms before this time, 
otherwise it would contradict local existence stated in the Appendix. □ 

The case of V = §(0, R) 

In the case where the set of admissible velocities is the sphere of radius R, we can adapt the proof above to demonstrate that 
blow-up occurs if 

S = R 2 M fxoRM_A (3. 7) 
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is a positive quantity. The important modification arises in the constants evaluated in Lemma fTSI We can adapt the computations 
in that Lemma to obtain that for any (p, q) 6 R 2 x K 2 , we have 

7ri? 3 

(p • v)(v ■ q)+ dv = —^-(P ■ q) ■ 
As a consequence, we get that the key differential virial inequality (|3.6|) becomes in this new setting, 

(x ■ v)f (x, v) dvdx ~St + xo2R 2 K{0) , 



(3.8) 



: xV 



This concludes the adaptation to the case of the sphere. 

3.2 Blow-up including chemical degradation (a > 0) 

It is now classical that for results concerning blow-up in the parabolic Keller-Segel system, the additional contribution of chemical 
degradation 

— AS + aS = p , 

does not change dramatically the flavor of the results. It affects the blow-up criterion, however, and so it does in our situation. 

In order to study the influence of the chemical degradation, we shall estimate carefully the size of the corrective terms that 
come from the difference between the Poisson kernel and the Bessel kernel, respectively 



1 



1 



B (z) = — log — , B a (z) = — 



+ 00 



1 



47T./ t=0 t 



dt. 



2tt "° \z\ 

Lemma 15. There exists a universal constant C such that 

||V-B a -V-B ||oo <V^C. 

Remark 16. Observe that the scaling ^/a comes naturally from knowing the estimate between B\ and Bq. 
Proof. We give an argument in Fourier space. In fact B a (^) = (a + |^| 2 ) -1 , and we get: 

1 

2^ 
2^ 
2^ 

\f0L 



(3.9) 



(3.10) 



\VB a - VBol 



< 



< 



< 



WB a - WB 

a 



< 



2tt 



lel(« + ie| 2 ) 
1 



1 



di 



ICI(i + ICI : 



dC 



□ 



Proof of Corollary Following (|3.2|) and subsequent lines we are able to obtain a perturbated virial identity for the second 
space moment of the density: 



d 2 T , N 
dP m 

ttR 4 



\v\ 2 f(t,x,v)dxdv + xo 



nR 4 



2R 3 

x ■ V5(i, x)p(t, x) dx — xq— ^— / / {x ■ f)|V S\(t, x)f(t, x,v) dvdx 



+Xo^M x -(vS(t,x)-VS(t,x))p(t,x)dx-xo— / / (x ■ v) (| VS\(t, x) - \VS\ (t, xfj f(t, x, v) 



(3.11) 



) dvdx , 



where X7S is defined by — AS — p as above (cf. the case a = 0). The first line of (|3.11[) can be explicitly computed as before. It 
remains to estimate the error terms (on the second line) using Lemma [T5l In fact we have using Young's inequality: 



Xo^- / x-(yS{t,x)-VS(t,x))p(t,x) 



dx 



< XoM- 



7rR 



\x\\\VB a -VB \\ oo p(t,x)dx 



kR 



< s/olxqM——C I \x\p{x)dx 



^R 



< ^x M^CVM v / 2J(t). 
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In the same way we get, 
2R 3 



Xo— / (x-v)(\VS\(t,x)-\VS\(t,xj)f(t,x,v)dvdx 



2R A 

< Xo— I I ' 

2i? 4 

< ^/axoM—^-C I \x\p(x)dx 
r 2R 4 



VS(t,x) -VS(t,x) p(t,x)dx 



Therefore we end up with the following intcgro-differcntial inequality instead of (|3.6p : 
c i r r 2R 3 

T I(t)< (x-v)fo(x,v)dvdx-5t + xo^ r K(0) + ^XoM 3 / 2 R i C y^dr. 
Recall that 5 > is defined in (13. 4| . The universal constant C is now fixed for the rest of this proof. Denote 



(3.12) 



Po 



2R 3 

(x ■ v)f Q (x, v) dvdx + Xo^-K(O) , 



V = V^XoM 3/2 R 4 C . 
Integrating once more and inverting the order of integration we get 

J(t) < 1(0) + p t -f 2 +V 



drds 



s=0 Jt=0 
t 



< 1(0) + p t - -t 2 + Tj {t-T)y/l{j)dT 

* Jt=0 

< 1(0) + p t - + vj (|(* - rf + 1/( T ) j <h 



S .2 , ^,3 , V 



< 1(0) + p t - -t l + -^t 3 + 



We obtain as a consequence, 



f I(t) dr < e^' 2e f e- ,1T/2e (l(0) + p 

t=0 Jt=0 V 



A sufficient condition for a contradiction to occur is that the right-hand side becomes negative as t — > oo. For this purpose, 
compute 



2emo 2e 2 <S 2 4e 4 

/(0) + CT — (7 + — 

77 jyZ j^Z 



2 1 a 3 ) da 



Choose for instance e = VS/2. Then the quantity J* =Q I(T~)dT eventually vanishes if the following condition is fulfilled: 

J(0) + ^< * 
r\ 2r\ l 

To conclude the proof of Corollary [31 observe that, resulting from (|3.5[) we get: 

,„<^y27(0) + ^ 3 



(3.13) 



so that the necessary condition (|3.13| might be replaced by the stronger criterion 

,2/(0) 



if 



2 RVM 



X oM 3 / 2 R z 



tin 
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This criterion can be read as X + 2A\[X < S, which is equivalent to \[X < V 6 + A 2 — A. Therefore, the following criterion is a 
necessary condition for solutions to blow-up after finite time, 



,2/(0) „ A2 S , o / , S 



a = rSm + *° m3/2r3 



3tt 

□ 

3.3 Drift-diffusion limit and blow-up (case a = 0) 

In the Introduction, we derive formally the parabolic Keller-Segel system from the kinetic system with suitable turning kernel 
Repeating the virial computation in this case leads to 



e 172 J ^(*) = / \v\ fe(t,x,v)dvdx + - (x-v)p e {t,x)F{v)dvdx-- x-j £ (t,x)dx-^—M z -xo^ 1 rK e (t) 
dt 1 J„ e ,L J„ e .L, 32 3 df 



> X J V 

?4 



because J v vF(v) dv = by assumption. To conclude as above that the remaining term is the derivative of a nonpositive quantity, 
observe that 



(V\x\ 2 ) Je (t,x)dx = / |x| 2 V- j e (t,x)dx 
Jm 2 

= -e \x\ 2 d t p e (t,x)dx 



d i 2 



-€— I \x\ A p 6 (t,x)dx. (3.14) 



Integrating once in time we obtain 



2 d r ^ ^ /* .Air /_ .A , n2w ft XoR 2 M\ t T ^ 2R 3 



ej t Ie(t)<e j J {x-v)f t {x,v)dvdx + R z M\l- ^\ i2 j t + Z e (0) - I e (t) + Xo-jj- (K e (0) - K e {t)) ■ (3.15) 

Arguing as before, we conclude that the solution blows- up in finite time under the same assumption as Theorem [2l 

On the other hand, we know precisely the blow-up criterion for the parabolic limit. It depends upon the choice of the 
relaxation function F(v). This function can vary between two extremal choices: a Dirac mass at zero, and a Dirac mass spread 

on the sphere = R}. Consider for example F(v) = j^-ly- The parabolic limit (|1.20|) writes 

d tP = v • v„) - • ( P VS) , 

thanks to isotropy (|3 . 1[) . Thus the blow-up criterion for finite-time blow-up reads 

1G 



M > 



xoR 2 



which differs from the kinetic criterion (II. 4|) by a factor 2 (the two criteria actually match in the case where F(v) is the normalized 
Dirac mass on the sphere of radius R). To fill this gap it is necessary to reconsider the parabolic limit a e — > 0. In fact the 
positive contribution R 2 M in (|3.15j) comes from the upper-bound of ff R2x y \v\ 2 f e (tx,v)dvdx. One may notice that for small e, 
f e (t,x,v) gets close to p(t,x)F(v), and this upper-approximation is not sharp (except in the case where F(v) is especially the 
Dirac mass on the sphere of radius R) . Therefore we may replace Jj R2 xV \v\ 2 f e (t,x 1 v)dvdx hy M J v \v\ F(v) dv. To finish with, 
observe that the diffusion tensor at the parabolic limit can be calculated due to rotational invariance, 



v ® vF(v) dv 



v 



~ U\v\ 2 F(v)dv^ Id. 



Under these considerations, the kinetic and the parabolic criterions do coincide. 
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4 Global existence for small mass 



The aim of this section is to prove global existence for the kinetic model (|1.3[) in the spherically symmetric case under the small 
mass condition stated in Theorem [4] 

Spherical symmetry is used in a crucial way in the following estimate: 





< 















M 
2tt\x\ 



(4.1) 



To justify (|4. 1[) . just write 



-S'(r) = 



Xp(X) dX < 



A=0 



M 
2nr 



Recall that < 7 < 1 is an exponent given by some upper-bound on the initial data. Introduce the auxiliary function 



k(x, v) = k{) 



= h) 



\n v ^(x)\- 



if (v ■ x) < 
if (vx)>0 



(4.2) 



where U v ± = Id — v ® w/|d| 2 denotes the orthogonal projection onto v ± . We shall prove in fact that, as soon as fo(x, v) < k(x, v), 
then it holds true that f(t,x,v) < k{x,v) for all time t > 0. This will be achieved through a comparison principle adapted to 
our context. 

Proposition 17 (Properties of the auxiliary function k). (i) The function k(x,v) belongs to L^ ocx L1 provided that p < 2/7 
and q < I/7. Moreover it is a C 1 function of (x,v) in (M 2 x V) \ {{x,v)\x TT v). 

(ii) For all (x, v) we have k(x, v) > fcol^l -7 - Therefore the initial comparison fo(x, v) < k(x, v) is guaranteed by the assumptions 
of Theorem^ 

Proof. Working exactly as in the proof of (|4.4[) below we find that J v k(x, v) q dv — fc^^-f^g) \x\ _79 , where fl(jq) (see (|1.6| l) is 
finite thanks to q < I/7. Taking the L^-norm gives the first assertion. The inequality k(x, v) > ko\x\ ~ 7 is obvious when x-v < 0, 



and follows from |Hy± (x)\ = \x\ \ sin6\, where 9 = Z(x, v), when x ■ v < 0. 

The crucial Lemma, which motivates the definition (|4.2| of k(x, v) is the following one. 
Lemma 18. Assume (jl.7[) . The function k(x,v) is a supersolution of (|1.3p . in the sense that: 



v ■ W x k 



V ■ 1 — r I -, — rk(x, v) 



1 M 



^ Xo[v- — ) -tt— I k(x,v')dv 



Proof. First we evaluate 



□ 



(4.3) 



k(x, v ) dv' — ko\x\ 7 



{v'ev\(x-v')<o} 



k \x\-^ + k \x\-< 



dv + fcoM - 

R ptt/2 

v =o Je=-iT/2 



f 


Il v ,x (x) 


>{v'£V\(x-v')>0} 


\x\ 



dv' 



sin 6\ vdvdQ 



=fc |*rMfi(7) 

IVI 

where we have used part (ii) of Proposition [171 

In order to prove (14. 3j) . let us distinguish between (v ■ x) > and (v ■ x) < 0. In the former case we have 



v ■ \7 x k 



-jk v ■ (IL V ± o U v ±) (x)\U v ± (x) 



-7-2 



0. 



(4.4) 
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because H v ± is a linear symmetric operator whose image is orthogonal to v. In the latter case (v ■ x) < we have 

v-\7 x k — — 7fco(w • x)\x\ 1 ~ 2 



7 



' ' " r I -r-Mx,v) 
x\ / \x\ 



> v 



2 7 1 



k(x, v') dv . 



\x\J_ |V|fi( 7 ) |ar| 

and we conclude the proof by using the smallness condition (jl.7p in the assumptions of Theorem |4j 
Definition 19 (Set of admissible exponents). A couple of exponents (j>, q) is said to be admissible if it satisfies 

2 1 „ 1 1 1 <f 1-7/ 2 

2<p<-, Kq<-, 0< <-, ^>— ' 

7 7 q p 2 p ' 1 



□ 



7 

T/iis sei is nonempty as it can be seen when (p, q) — > (2 + , 1" 



' 7 



Lemma 20 (l^.L q regularity is ensured by comparison). Let (p,q) be a set of admissible exponents. Assume that f(T,x,v) lies 
below k(x, v). Then f(T, x — tv, v) belongs to L^L q for all t > 0. 



Proof. We shall first prove that 



For this purpose, let us decompose 



Vx^O \\k(x-tv,v)\\ Ll <C(q,^k ,V)\x\ 



(4.5) 



k(x-tv,v) q dvk q / \x- tv\~~> q dv + k q / \U V ± (x - tv)\^ q dv . (4.6) 

V J {v£V\(x-tv)-v<Q} J {v€V\(x-tv)-v>0} 

For the second contribution in the right-hand-side above we use Tl v ±(x — tv) = 1X^(2) and |n„^ (a:) | = \x\ |sin#|, where 
9 = Z(x, v), to get 



k q 



{veV\(x-tv)-v>0} 



\Ii v ±{x - tv)\-^ q dv < k q I \U v ±(x)\^ q dv 



2 Je=—K 
= C(q, 1 ,k ,V)\x\- iq . 



(4.7) 



To estimate the first contribution in the right-hand-side of (|4.6|) we distinguish between two cases. 
If \x\ < 2tR then V C {v e JR 2 | \x/t - v\ < 3R}, therefore 



\x-tv\-~< q dv < k q t-^ q 
{veV\(x-tv)-v<0} -i \ 



-79 



dv 



< kit'"" 1 



{tu€R 2 |M<3-R} 



< k q 



2R 



-~<i 



C(q,~f,V) 



< C(q, 7 ,k ,V)\x\-^ . 
On the other hand, if |x| > 2tR then we use \x — tv\ > |ar| — t \v\ > \x\/2 to obtain eventually, 



k q 



\x-tv\-~< q dv < k q 2 79 |x|" 79 \V\ 



(4.8) 



' {veV\(x-tv)-v<0} 

In a second step we split the T v x L q norm of f(T, x — tv, v) into a short-range part (in space) and a long-range part, as follows 
\\f(T,x-tv,v)\\ L p L « < \\f(T,x - tv,v)l\ x \<i\\ L p Li + \\f(T,x - tv,v)l\ x \>i\\ L p L , . 
For the short-range contribution \x\ < 1 we use (|4.5j) - which is a combination of (14. 7|) and (14. 8|) - to get 



|/(T,x- tv,v)l\ x \ 



< \\k(x- to,w)l| a .|< 1 || iPi9 < C(q,j,k ,V) \x\ 7 1 



<C(p,g,7,Ao,V). 
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because jp < 2. 

For the long-range contribution \x\ > 1 we introduce a pair of auxiliary exponents (P, Q) such that, 

2 ^ 1 P' v ' 

- <P<oo, 1 < Q < - , 7^ = -- 

7 i Q q' 

We shall ensure that such a choice of (P, Q) exists: in fact when Q — ► (1/7)" we have 

1 - 7/2 < <i_ = , 1 1 



1-7 p' P' P' 1 - 7 
Therefore we can find P' < (1 - 7/2) -1 , i.e. P > 2/7. 

Define 8 = P'/p' = Q'/q'. We have the interpolation relation L%.L% — [L\ v , L^L^]i- 



because 



1 

- = 1 

P 



Q 



As a consequence, 



\f(T,X- *W,u)l|x|>l||jr,g £ g < to >' y ) ]L bl>l|LiLi ||/( T ' X_ to ' W ) 1 kl>l|lz J: F>i« 



< C'(p,<Z,7,fco,V,M) bpl^ 

< C(p,g,7,fco,V,Af). 



□ 



Proof of Theorem ^ We obtain from Lemma [18] and the elliptic bound (|4.1I) (which holds true in the spherically symmetric 
framework) the following crucial estimate: 



v ■ S7 x k > xo [v 



k(x,v')dv' > xo(v ■ VS)+ / k{x,v')dv' 



M 

\x\J _ 2tt\x\ J v , 

Therefore we get the following differential inequality, well-suited for proving a sort of comparison principle: 

d t (f-k)+v V,(/ -k)< X o(v ■ V5)+ (p(x) - / k(x, v 1 ) dv^ 



(4.9) 



(4.10) 



However the non-local nature of the right-hand side requires more regularity, and a local in time estimate (obtained in the 
Appendix) will enter into the game. Due to the lack of integrability at infinity of the reference function k, we aim to localize 

in space such a partial differential inequality. In order to do so, multiply (|4.10[) by the test function ip(x) = exp {^—y/l 

which satisfies: 



Introduce the notation: 



Vz \Vip(x)\ = 



K{x) 



VT+W 

k(x, v') dv' 



<p(x) < <p(x) 



(4.11) 



v'ev 



We get the following localized partial differential inequality, 

dt ((/ - k)<p) + v ■ V x ((/ - k)<p) < X o(v • V5)+ ( P -K)<p + (f- k)v ■ . 

Introduce P e a sub-approximation of the positive part (■)+. Namely we choose P e (z) — eP\{z/e), where 

if z < 

P^z) = ^ z 2 /2 if < z < 1 . 
z-1/2 if z>l 



(4.12) 
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In particular, P £ (z) is identically zero on {z < 0}, < ~P' e {z) < 1 and moreover 

Vz \z\P' e {z) < 2P e (» . (4.13) 

Multiplying (|4"7T2|) by P £ ((/ - we obtain 

d t P e ((f - k)<p) + v ■ V x P e ((/ - k)p) < X o(v ■ V5)+ (p - K) pP' E ((f - k)<p) + (/ -*)(«• Vp) P^((/ - k)<p) 

< X o(v.VS)+{p-K) + <p + 2KP e {(f-k)<p), 

due to (|4.1ip and (|4.13[) . Therefore we obtain, thanks to Duhamel's representation of the solution after a given time T, 

P £ ((/ - k)<p)(T + t, x, v) < e 2Rt P e ((f - k)<p)(T, x - tv, v) 



e 



mt X oR I e~ 2Rs (|V5| (p - K) + ip) (T + s,x — (t — s)v) ds . 



s=0 



As £ — + we obtain the following integral inequality, 

(ftp - k<p)+{T + t, x, v) < e 2Rt {f<p - k<p) + (T, x - tv, v) 



e 



2Rt X0 R / e- 2Rs (|VS| (pep - K<p) + ) (T + s,x — (t — s)v) ds . 



s=0 



Next compute similarly as in the Appendix (by the help of the dispersion estimate Lemma [6]) for admissible exponents (p, q) 
(Definition [111), 



„2Rt I 



+ Xo C(V)e 2Rt [ e-^"- r7 l_ 7 _|||V,S'|l/ V ■ K^) + iT ^ ,)\\ L ,,I, 

Js=0 

< e 2Rt \\(fp~k<p)+(T,c 

+ X0 C(V)e 2Rt [ e- ll: ~- ,;, .J IV.s-tr + Mll^, , f ,s - Ks) . , T - .s||| L , ,1, 

Js=0 



- tV,V)\\ L P L * 




-2Rs 1 




( t - s)2(l/9- 


-i/p) 


- tv,v)\\ L p Li 




~2Rs 1 




(t_a)a(x/9- 


-i/p) 


'- tv -, v ))\\LlLl 




-2Rs 1 




(t_a)a(x/9- 


-i/p) 



< e 2ffl ll(/ 

+XoC(F)e 2 * £ q e - 2Rs j r —^ ZW) \\p{T + sJMKw - A» + (T + s)\\ LP ds , 

where the Holder exponents are given by 1/p* + 1/p = 1/q and the Sobolev exponent (r < 2) is given by 1/r = 1/2 + l/p* = 
1/2 + l/« - 1/p < 1. 

Furthermore, because the positive part (•)+ is a convex function, and V is a bounded set, we have by Jensen's inequality, 

(p<p-K<p) + (T + s,x) = ( [ \V\(f<p-k<p)(T + s,x,v)^) <[ QV\(f<p-k<p)) + {T + s,x,v)^, 

\Jvev I v I / + Jvev \ v \ 

|| - A»+(T + ,s)|| LP < HC/V - + s)|| igi i < C(q, V)\\(f<p - ky) + {T + s)|| iM . 

As soon as /(T, x, v) remains below k(x, v), Lemma [201 guarantees that the free transport contribution f(T, x — tv, v) belongs 
to LI for admissible exponents (p, q). As a consequence of the local in time existence result of Proposition[22] (in the Appendix), 
we can ensure that f(T + t, x, v) belongs to L V X L1 for small time t > 0. Thus p(T + t, x) belongs to L r for any r < 2 in particular. 

Therefore the Gronwall lemma guarantees that ||(/y — kp) + \\ L P L 9, if it is zero up to some time T, it remains zero for small 
later times T + t. 

We have proven a comparison principle which prevents solutions to (| 1 . 3|) from blowing-up. □ 
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Normale Superieure for their hospitality and financial support during his sabbatical leave in the spring semester of 2008. 
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Appendix A: Solutions having spherical symmetry 

The notion of spherical symmetry in kinetic theory is contained in the following Definition 1211 Recall that the set of admissible 
velocities is the ball V = 6(0, R). 

Definition 21. A function f(x,v) defined for (x,v) £ R 2 x V is spherically symmetric if for every rotation of R 2 we have 
f{Qx,Qv) = f(x,v). 

If f(x, v) is spherically symmetric then the space density p(x) — J veV f(x, v)dv is spherically symmetric in the usual sense, 
i.e p(Qx) = p(x) for all rotations 0. Therefore p depends only on r = \x\. Abusing notations we write both p(x) and p(r) but 
the meaning will always be clear from the context. 

If f(x,v) is spherically symmetric then its current j(x) — J v vf{x, v)dv does not necessarily point in the direction of x. For 
example, if f(x, v) = x ■ v = x 2 vi — Xiv 2 , then, 

j(x) = -^-(-x 2 ,xi) ± x. 



However, if we decompose 

j(x)=j\x)-^ J +j ± (x)^, (4.14) 



x 

then we have for every rotation 0, 



j(Qx)= / vf(Qx,v)dv= / (Gw)f(Qx,Qw)d(Gw)Qj(x). 
Jv Jv 

Therefore the decomposition's coefficients in (|4.14l) are both spherically symmetric: 

jHQx)=jW(x), and ^(Qx) = j^{x) . 

Abusing notation we write 

j(x) = j"(r)- +j ± (r)— . 

r r 

We can then simply derive the following identities which will be crucially used in Section [3] 

x- j(x) =ri ll (r), (4.15) 
V-j(x) = ^(rjW(r)y . (4.16) 

The kinetic system (11.31) preserves spherical symmetry 

If we start with spherically symmetric initial data fo(x, v) then Proposition ! 1 01 guarantees the existence of a local in time solution 
f(t,x,v). It is not difficult to verify that for any rotation the function f(t,Ox,®v) is also a solution to (|1.3p (see below). 
Therefore, by the uniqueness part of Proposition 1101 thev have to coincide. It follows that f(t,x,v) is spherically symmetric 
throughout the time interval of existence. 

Let (/, S) be a solution and be a rotation of R 2 . Define (g, Q) by 

g(t, x, v) = f(t, Qx, Qv) , Q(t, x) = S(t, Ox) . 

On the one hand, 

d t g(t,x,v) + v-V x g(t,x,v) = d t f(t,@x,Qv) + v Q T (V x f)(t,Qx,Qv) 

= d t f(t, Qx, Qv) + Qv ■ (V x f){t, Qx, Qv) . 
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On the other hand, 

T[S](t, Qx, Qv, v')f(t, Qx, v') dv' - [ T[S] (i, Qx,v', Qv)f(t, Qx, Qv) dv' 

Jv 

= [ (Qv (VS)(t,Qx)),f(t,Qx,v')dv' - ( (v' ■ (VS)(t,Qx)) + f(t,Qx,Qv) dv' 

Jv 

(v ■ 8 T (V5)(i, 9a;)) f(t, Qx, Qw) dw ~ [ (w T (VS)(t, Qx)) f(t, Ox, Qv) dw 

Jv 

(v ■ (WQ)(t,x)) + g(t,x,w)dw - / (w ■ (VQ)(f, Qx)) + g(t, x, v) dw 

Jv 

T[Q](t,x,v,w)g(t,x,w)dw — / T[Q](t, Qx, w, v)g(t, x, v) dw. 
jv Jv 

Also 

-AQ(x) +aQ(x) = -(AS)(Qx) +aS{9x) = / f(Qx,v)dv= [ f{Qx,Qw)dw= [ g{x,w)dw 



Appendix B: Existence and uniqueness in a weaker framework 

The goal of this appendix is to provide a variant of Proposition [TO] in a framework well-suited for proving the global existence 
result of Section |4j As a matter of fact, the reference function k(x, v) used there belongs to Lf ocx Ll for any 1 < p < 2/7 and 
1 < q < I/7. On the other hand, Proposition [TO] deals with solutions lying in L V X L V V , p > 2, thus it does not cover properly the 
case 7 > 1/2. 

Assumption A2 (Initial datum, precised version). Assume that the initial density fo(x,v) > belongs to L\ v and satisfies the 
estimate ||/o(a; — tVyV)^^^ < Co for small times t>0, and for some couple of exponents (p,q) verifying: 

2<p, Kq, 0<i--<~. (4.17) 
q p 2 

Observe that Assumption A^]is satisfied as soon as /o € L v x v for some p > 2. Thus it is weaker than Assumption A[T] except 
for the condition of spherical symmetry. 

Proposition 22. Assume that the initial density fo verifies Assumption A\^for some couple of exponents {p, q) such that (|4. 17j) 
holds true. Then there exists a unique (local in time) solution to system (|1.3p with f(t,x,v) £ I^.h\. 

Before the proof of this Proposition, let start with a general feature of kinetic transport equations with source and decay 
terms. The solution of d t h + v ■ V x h + \(t, x, v)h = g with vanishing initial data is given by the Duhamel's representation, 

h(t,x,v) — J g(s, x — (t — s)v, v) exp |— J A(r, x — (t — r)v, v) drX ds . 

In case A is a nonnegative function, we obtain, 

\h(t,x,v)\ < J \g(s, x — (t — s)v, v)\ exp /— J A(r, x — (t — t)v, v) rfr| ds 

< / \g(s,x- (t-s)v,v)\ ds. 

Js=0 

Proof. We aim to write directly a fixed-point argument under the reference norm, 

\\9\\y= sup (\\g(t,x,v)\\ „ +\\g(t,x,v)\\ L , Li ) . (4.18) 

0<t<T K ' 

We start from 

d t (fi -f 2 )+v- V a (/i - h) + Xolo |VS 2 | (/1 - f 2 ) = 

Xo ((v ■ V5i) + - (v ■ VS 2 ) + ) Pl + X o(v- V5 2 ) + ( Pl - p 2 ) - xo" (\VS 1 \ - \VS 2 \) fi ■ (4.19) 
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Applying the preliminary observation to equation (|4.19D we obtain 

I/1-/2IW < C J ^\\VS 1 -VS 2 \\ 0O \p 1 (s,x-(t-s)v)\ + \\VS 2 \\ 0O \p 1 -p 2 \(s,x-(t-s)v))ds 
+C f \\WS 1 -WS 2 \\ OQ \f 1 ( s ,x - (t - s)v,v)\ ds, 

Js=0 

Therefore we are able to develop a dispersion technique as usual, 

\\(fi-h)(t)U Li < c f (t - s )-2(i/ 3 -v^) ( ikvSx - v^x^iu + Hv^h^ikpx - P2 )( S )|U«) 



ds 



s=0 

t 



+C / ||(V5i-V5 2 )(«)|| 00 ||/ 1 (*,a;-(t-s)t;,t;)||iPi,da 

Js=0 

< CQ\fi\\ Y + \\M Y ) (^ Q (*- S )" 2(1/ ^ 1/P)rfs ) Wh-hh 

+c( f II fx (s, x -(t-s)v,v) W&L'ds) \\fi-f 2 \\Y, 



\Js=0 J 

where we have used Lemma ITT1 In parallel, we get a bound for ||/i(s, x — (t — s)v, v)\\ L p . We argue as follows: f\ solves 

dtfi + v ■ V x fi + xou \VSi\ ft = xo(v> V5i) + P i , 
and since the coefficient \qlo |VSi| is non-negative, we can argue as above to get 

\fi(s,x,v)\ < \h(0,x- S v,v)\+C f WVS^W^lpii^x - (s - t)v,v)\ dr, 

Jt=0 

and eventually 

\fx(s,x-(t-s)v,v)\ < |/i (0,x I +C f ||VS , i(r)|| 00 | / 9i(r,a;-(t-r)«,i;)| dr , 

Jt=0 

x-(t- s)v, v)\\ l%l p < \\f (x - tv,v)\\ L P Li +C f (t - rr^-^llVMr)!!^ \\ Pl (r)\\ L , dr 

Jt=0 

<C + C||/i||£ / (i-Tj-^'-Vrtdr, 



from which we deduce that the flow is contractant for small time w.r.t. the reference norm || • \\y, and relatively to the initial 
datum. 

□ 
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